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Abstract 

In this paper, we show how concentration inequalities for Gaussian quadratic form can 
be used to propose exact confidence intervals of the Hurst index parametrizing a fractional 
Brownian motion. Both cases where the scaling parameter of the fractional Brownian motion 
is known or unknown are investigated. These intervals are obtained by observing a single 
discretized sample path of a fractional Brownian motion and without any assumption on 
the parameter H. 
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1 Introduction 



Since the pioneer work of Mandelbrot and Nessl ( 19681) . the fractional Brownian motion (fBm) 
has become widely popular as well as in a theoretical context as in applications. Fractional 
Brownian motion can be defined as the only centered Gaussian process, denoted by (Bjj(t))teB., 
with stationary increments and with variance function «(•), given by v(t) = C 2 \t\ 2H for all 
t £ 1. The parameter H € (0, 1) (resp. C > 0) is referred to as the Hurst parameter (resp. 
the scaling coefficient). In particular, when H = 1/2, it is the standard Brownian motion. 
In general, the fractional Brownian motion is an iJ-self-similar process, that is for all 5 > 0, 
(Bh (<5i)) tgR = 8 H (Bu(t)) t ^ (where = means equal in finite-dimensional distributions) with 
autocovariance function behaving like 0(|fc| 2 ) as \k\ +oo. Thus, the discretized increments 
of the fractional Brownian motion (called the fractional Gaussian noise) constitute a short- 
range dependent process, when H < 1/2, and a long-range dependent process, when H > 1/2. 
The index H characterizes also the path regularity since the fractal dimension of the fractional 
Brownian motion is equal to D = 2 — H. G ene ral references on self-s imilar processes and long- 
memory processes are given in lBeran f|l994h or lDoukhan et al.l (|2003h . 

The aim of this paper is to propose confidence intervals for the Hurst parameter based on 
a single observation of a discretized sample path of the interval [0, 1] of a fractional Brownian 
motion. To do so, the most popular str ategy consists in u sing the asymptotic normality of some 
estimators of th e Hurst paramete r, see Coeuriollvl (120001) for a survey on the estimation of the 
self-similarity or Shen et al. ( 2007 ) and Coeuriollvl ( 20081 ) for more recent discussions in a robust 
context. R ecently, a new strat e gy ba sed on concentration i nequalities for Gaus sian processes 
obtained bv iNourdin and Viens (|2009l) has been proposed bv iBreton et al. (2009). In this case, 



the confidence intervals are non- asymptotic and they appear to be very interesting when the 
sample size is moderate. Our contribution is to improve this direction both from a theoretical 
and practical point of view. In order to present our different contributions, let us first recall the 
confidence interval proposed bv lBreton et al. I (120091) . 



Proposition 1 Assume that one observes a fractional Brownian motion at times i/n for i — 
0, . . . , 71+ 1 with scaling coefficient C = 1 and with Hurst parameter satisfying H < H* for some 
known H* 6 (0, 1). Fix a 6 (0, 1), then for all n large enough satisfying q n (a) < (4 — A H )\/n, 

^log(|), 



where q n {a) := \ (b(a) + \Jb(a) 2 + 852 log (§)J with 6(a) : 



have 



H E 



> 1 



(1) 



where for t > 



9n (K nf (t)) 



log 1 - 



(4-4«*)V^ 



1 log(S„ 

2 21og(n) ^ 21og(rc) 

1 logQgn) lQ g( 1+ (4-4"*)^) 

2 2 log(n) 2 log(n) 



where g n is the function defined by g n (x) — x — ^°2\o S (n) ^ ant ^ is the following statistic 



(2) 
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Let us give some general comments on this result. First, note that this procedure cannot 
be applied to a fractional Brownian motion whose scaling coefficient C is unknown. Secondly, 
important drawbacks of this procedure rely upon the assumptions made on H* and n, which 
exclude the possibility to use this confidence interval when the sample size is small: 

• Given a and H*, the following table presents the minimal value of the sample size n in 
order to ensure that q n (a) < (4 — 4 H )y/n. 







0.1 


0.2 


0.3 


0.4 


H* 
0.5 


0.6 


0.7 


0.8 


0.9 


a 


= 1% 


271 


298 


335 


388 


471 


611 


886 


1592 


4936 


a 


= 5% 


189 


208 


233 


270 


328 


425 


617 


1108 


3437 


a 


= 10% 


154 


169 


190 


220 


266 


346 


501 


900 


2791 



• The following table exhibits the maximal value of H* , denoted by H* , required in order to 

ensure q n (a>) < (4— 4 H )y / n in terms of a and n. Note that H* = log (max (1, 4 — q n (a)/y/ri)) / log(4), 
which means that, given a and n, a confidence interval is only available for H G (0, H*). 











n 












50 


100 


200 


500 


10000 


10000 


a 


= 1% 


0.00 


0.00 


0.00 


0.53 


0.93 


0.93 


a 


-5% 


0.00 


0.00 


0.17 


0.65 


0.94 


0.94 


a 


- 10% 


0.00 


0.00 


0.34 


0.70 


0.95 


0.95 



We are now in position to specify our different contributions: 

• We sli ghtly improve the bounds of the concentration inequality obtained bv lNourdin and Viensl 
(|2009h . see Section [2] and Proposi tion [2] for more details. Note in particular that, in con- 



trast to lNourdin and Viens (|200£ 
optimize numerically our bounds. 



and 



Breton et al ] (l2009h . we are tracing the constant to 



In the case where the scaling parameter C is known, we propose a new confi dence interval 



witho ut any preliminary assumption on the Hurst parameter H (in contrast to lBreton et al 



(120091) 1 and with a very slight condition on the sample size. For instance, in comparison to 
the previous tables, our confidence interva l is computable as soon as n > 3. Furthermore, 
by using ideas similar in ICoeuriollv ( 2001 ) for the problem of the estimation of the Hurst 
parameter, we also propose a confidence interval when the scaling parameter C is unknown. 
This new confidence interval has the nice property to be independent of C and independent 
of the discretization step. It is remarkable that, in the both cases (C known or unknown), 
the lengths of the confidence intervals we propose behave asymptotically like the ones 
derived in an asymptotic approach, that is they behave like l/y / nlog(n) when C is known 
and 1 / Jn when C is unknown. 



As suggested by the expression of the statistic in (0) , the procedure described in Proposi- 
tion Q] is based on the increments of order 2 of the discretized sample path of the fractional 
Brownian motion. Taking the increments of order 2 is a special case of filter to work with 
and it is known that discrete fi ltering has been propo sed a nd used in an est imation con- 
text, see llstas and Land (|l997t ). lKent and Woodl (|l997t ) and lCoeuriollvl (|200lh . Recall that 
the main interest in filtering the fractional Browian motion is that the action of filtering 
changes the correlation so that, for instance, the increments of order 2 of the fractional 



3 



Brownian motion constitute a short-range dependent process (i.e. its correlation function 
is absolutely summable). Such a behaviour is required to obtain an efficient concentration 
inequality. In this paper, we propose to construct confidence intervals not only based on 
the increments of order 2 but on more general filters such as, for instance, increments of 
larger order or the Daubechies wavelet filters. . . Finally, let us also underline that a crucial 
step consists in obtaining an upper-bound of the supremum on the interval (0,1) of the 
i 1 — norm of the correlation function of the discrete filte red series of the fract ional Brownian 
motion. When considering the increments of order 2, iBreton et al ] (l2009h have obtained 



the bound 17.75/ (4 — 4 H ). We have widely improved this point since we compute explicitly 
this supremum for a large class of filters (including increments of order 2). As an example, 
for the increments of order 2, this gives the explicit value 8/3. 

• Based on a large simulation study, we assess the efficiency of the different procedures that 
we propose and we compare them with ones based on an asymptotic scheme. We discuss 
and comment these results. 



The rest of this paper is organized as follows. In Section^ we give the concentration inequal- 
ities specially designed for our purposes. The filtering setting is introduced in Section [3] where 
the bounds for the £ 1 -norm of the correlation function of the filtered series are also obtained. 
Our confidence intervals for the Hurst parameter are proposed and proved in Section 2J both 
when the scaling parameter is known or unknown. Our results are discussed and compared to 
the literature in Section [5] Finally, computations expliciting some bounds for some special filters 
are given in Appendix [21 



2 Concentration inequalities 



Proposition[T]above is based on concentration inequalities proposed bv lNourdin and Viend (l2009h 



(see Proposition [3]) for smooth enough random variables with respect to Malliavin calculus (see 
Theorem 4.1-i)). By applying such inequalities to the random variables y/nV n where V n = 
i Y^7=i Hzi-X-i), H2(t) = t 2 — 1 is the second Hermite polynomial, and X = {-Xj}i<,-< n is a 
stationary Gaussian process with variance 1 and correlation function p, we obtain concentration 
inequalities for i^— variations of stationary Gaussian processes. In the sequel, for a sequence 
(ui) i6Z) we set := Y,\i\< n N- 

Proposition 2 Let K n — 2\\p\\(i . Then, for all t > 0, we have: 

P (VnV n >t) < <p r , n {*> K ") := f 1 + ( 3 ) 

P(V^V n <-t) < w,n(*;«n):=e^ (l--j=Y n l [(W Sl(<)- (4) 

Note that Proposition [2] can be applied to short-memory as well as to long- memory stationary 
Gaussian processes (as soon as n remains finite) . In order to derive Proposition [5] below, we 
shall briefly use some no tions of Malliavin c alculus. We just recall the only necessary for our 
argument and we refer to IBreton et al. I (j2009h and references therein for any further details. We 



stress that, once Proposition [5] is derived, only basic probability tools will be used. Without 
restriction, we assume the Gaussian random variables Xi have the form Xi = X(hi) where 
JT(K) = {X(h) : ft, e H} is an isonormal Gaussian process over a real separable Hilbert space H 
and {hi : i — 0, . . . , n} is a finite subset of H verifying ¥.[X(hi)X(hj)] = p(i—j) = (hi, hj)n. With 
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such a representation, V n can be seen as a double Wiener-Ito integral with respect to X, i.e. 
V n = h yi"-n hj g> hj) . In the sequel, to make easier the presentation, we rewrite Th. 4.1 of 
Nourdin and Viensl ( 20091 ) only for such random variables, see Proposition [21 Actually, in order 



to optimize our forthcoming results, Proposition [3] is a slight improvement of Th. 4.1. Before, 
recall that multiple Wiener- Ito integrals I q {f) are well defined for / e K 09 , the gth symmetric 
tensor product of H, q € N\ {0}; the Malliavin derivatives D transforms random variables (in its 
domain) into random elements with values in H; multiple Wiener-Ito integrals are in the domain 
of D and we have D t (I q (h)) — ql q _i(h(-,t)). Recall also that the Hermite polynomials H q are 
related to multiple Wiener-Ito integrals by H q (Ii(h)) = I q (h® q ) when \\h\\x = 1; in particular, 
for q = 2, we obtain h(hf - 1 = I 2 {h® 2 ). 

Proposition 3 Let Z = I%{f) satisfying 

\\DZ\\l < aZ + b (5) 

for some constants a > and b > 0. Then, for all t > 

at 

T 



\Z >t) < tp r (t;a,b) := e~~ 1 



\Z<-t) < ipi(t;a,b):=e- [1--) l [0 , 6/o] (t) 



(it 



Proof: The proof is a slight improvement of the bounds in (jNourdin and Viensl . I2009L Theorem 
4.1) obtained by a careful reading of the proof (with the following correspondance with the 
notation ther ein: gz(Z) = hWDZ ] ]^, a = a/2 and (3 — b/2). Denoting by h the density of Z, the 



argumen t of ( Nourdin and Viensl. 120091 Theorem 4.1) is based on the following key formula (see 



(3.16) in lNourdin and Viens (2009)) 



For the sake of self-containess, we sketch the main steps of the argument. For any A > 0, 
define m A : [0, +00) -)• K by m A {6) = E [e 9Z l { z<A}] ■ We have m' A {e) = E [Ze 0Z l {z < A} ] and 
integration by part yields 

m' A (6) = / xe 9x h(x)dx 

J — OO 

< 9 f A e 6x ( yh(y)dy) dx (7) 



p+00 

A 

we obtain for any 9 € (0, 2/a): 



< -E W\DZ\\ie ez l {z < A} ] . (8) 

where ([7]) comes from J A °° yh(y)dy > since E[Z] = 0, and (JSJ) comes from ((B|). Because of ([S]), 

2/a): 

m' A (0) < -^- mA {6). (9) 
2 — 9a 

Solving ©, using m A (0) = ¥(Z < A) < 1 and applying Fatou's Lemma (A — s- +00) yield the 
following bound for the Laplace transform and any 9 € (0, 2/a): 



b „ 2b 



E\e ez ] < exp ( -~6 - ^ In I 1 - 

a \ 2 



at 
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The Chebychev inequality together with a standard minimization entail: 

mm {- \t + - )6 -In 1 

K ee(o,2/a) { \ a J a 2 \ 2 



The minimization is achieved in 9 = (2t)/(at + b) and gives the first bound in Proposition [3] 
Applying the same argument to Y = —Z, satisfying ||£>V||^ < — aY + b, we derive similarly the 
second bound. Note in particular that condition [5] implies that Z > — b/a so that the left tail 
only makes sense for t G (—b/a, 0). □ 



Remark 1 Nourdin and Vienk \200& ) have obtained the bounds 



4>i(t; a, b) = exp — —J and 4> r (t; a, b) = exp 



TaWeQ] proposes a comparison of these bounds with ours through the comparisons of the values of 
their reciprocal functions since these quantities are of great interest for the considered problem. 
Observe that the most important differences occur when n is moderate. The example a = 4/i/n 
and 6 = 4 corresponds approximately to the choices of parameters that will be used in the next 
sections. 









a — 


1% 


a = 


2.5% 


a = 


5% 


a = 


10% 












<pT\<*) 




<pT\<*) 








11 = 


50 


NV 


6.0697 


9.2102 


5.4324 


7.9062 


4.8955 


6.8751 


4.2919 


5.7878 






BC 


4.4720 


7.1547 


4.1398 


6.9040 


3.8372 


6.0847 


3.4712 


5.2008 


n 


100 


NV 


6.0697 


8.1851 


5.4324 


7.1048 


4.8955 


6.2383 


4.2919 


5.3107 






BC 


4.9090 


7.3551 


4.4966 


6.4575 


4.1314 


5.7249 


3.7012 


4.9267 


ii = 


500 


NV 


6.0697 


6.9492 


5.4324 


6.1322 


4.8955 


5.4606 


4.2919 


4.7235 






BC 


5.5334 


6.6309 


5.0017 


5.8810 


4.5449 


5.2591 


4.0218 


4.5708 


n = 


1000 


NV 


6.0697 


6.6801 


5.4324 


5.9190 


4.8955 


5.2891 


4.2919 


4.5930 






BC 


5.6877 


6.4641 


5.1259 


5.7478 


4.6462 


5.1513 


4.1000 


4.4883 


n = 


10000 


NV 


6.0697 


6.2567 


5.4324 


5.5819 


4.8955 


5.0168 


4.2919 


4.3850 






BC 


5.9475 


6.1931 


5.3345 


5.5312 


4.8159 


4.9757 


4.2308 


4.3536 



Table 1 : Computations of the quantities <p l and (p r 1 (a) for the bounds obtained by 



Nourdin and Viensl (|2009f) (NV) and ours (BC) (see Remark [T] and Proposition [3]) for differ- 



ent values of n and a and for the particular case where a — 4/^y/n and b — 4. 



Remark 2 Note that tp r (-;a,b) (resp. tpi(-;a,b)) is a bijective function from (0,+co) (resp. 
(0,b/a)) to (0,1). Obviously, the index I in ipi (resp. r in ip r ) indicates we consider the left 
(resp. right) tails. 

We explain now how Proposition [5] derives from Proposition [3] standard Mall iavin calculus 



shows that, for Z = -JnV ni \DZ^ = — V"^ 1 X(i)X(i)p(i — i). see Theorem 2.1 in lBreton et al 
(|2009h . The following lemma ensures that condition (J5]) in Proposition [3] holds true with a 
2K n /\/n and b = 2n n . 



Lemma 4 For Z = \/nV n; we have \\DZ\W < n n (^j^Z + 1 

The proof of Lemm alHis a very slight modification of the first part of the proof of Theorem 3.1 in 
Breton et al. (2009) to which we refer. Finally, Proposition [3] applies and entails Proposition [3] 
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3 Applications to quadratic variations of fractional Brow- 
nian motion 



3.1 Notation 

From now on, Bh stands for a fBm with Hurst parameter H £ (0, 1) and with scaling coefficient 
C > and B# is the vector of observations at times i/n for i = 0, . . . , n— 1. We consider a filter 
a of length £+1 and order p, that is a vector with £+1 real components a,, < i < £, satisfying 



9=0 



q J a„ 



for j = 0, . . 



. ,p — 1 and (? p a 9 7^ 0. 
9 =o 



(10) 



For instance, we shall consider the following filters: Increments 1 (a = {—1, 1} with I = 1, p = 1), 
Increments 2 (a = {1,-2,1} with ^ = 2, p = 2), Daublets 4 (a = {-0.09150635, -0.15849365, 
0.59150635, -0.34150635} with i = 3, p = 2), Coiflets 6 (a = {-0.05142973, - 0.23892973 
0.60285946 -0.27214054, -0.05 142973, 0.01107027} with £ = 5, p = 2), see e.g. IPaubechiesI 
(2006) and iPercival and Waldenl (|2000h for more details. Let V a denote the vector B# filtered 
with a and given for i = £,..., n — 1 by 



r" 



9=0 



Let us denote by 7 r?y(-) and ptr(-) t he covariance and the correlation functions of the filtered 
series given by fsee lCoeuriollvl (120011) 1 



1 { 

E[V a (k)V a (k+j)] = C 2 XTt a H (j) with ^ = --^2 a q a r \q-r + j\ 



2H 



9,r=0 



and p%{-) ■= itfCO/^j^O) which is independent of C. Finally, define and V® as 



(11) 



^ n— 1 

S" 2 := V V a 



and 



,2ff 



71— 1 

w ■ C 2 ^(0) " 



,2i? 



c 2 ^(o) x ya U 



Note that V£ = S™=~/ ^2 ) where #2 (t) = t 2 — 1 is the second Hermite polynomial and 
X a is a stationary Gaussian process with variance 1 and with correlation function p a H . Observe 
that V®, n > 1, satisfy a law of large number (LLN) and a central limit theorem (CLT) 



(12) 



with explicit variance cr 2 Ja , see Proposition 1 in Coeuriolivl ( 2001 ). used to derive standard 
confidence interval for H . In contrast, our argument relies on concentration inequalities: applying 
Proposition [2] with these notation, we obtain fo all s,t > 0: 



(-s < Vn - £V° <*)>!- (p r , n -t(t; <,h) ~ <Pl,n-t(s; <_ ff ) 



(13) 
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where h = 2^, i , <n As previously explained, the action of filtering a discretized 

sample path of a fBm changes the correlations into summable correlations for the i ncrements. 
More precisely, it is proved that, for some explicit kn, Ph^) ~ \i\ 2H ~ 2p , see e.g. ICoeuriollv 



(120011 ). Thus, p a H {-) is summable if p > H + 1/2, i.e. p a H (-) is summable for all H G (0, 1) for 
p > 2 and only for iJ G (0, 1/2] if p = 1 (in the case i? = 1/2, observe that p^ 2 (k) — for all 

|*|>*). 

One of the aim is to obtain bounds in (|13[) independently of H and easily computable. Since 
fi,n(t, •) and tpr,n(t) ') are non-decreasing, the bound remains true with k" := 2 sup ffg ( ,r) ll/ 9 ?? Il^ 1 (z) 
replacing n n ,H- Here, and in the sequel, we set t = 1/2 when p = 1 and t = 1 when p >2. The 
following section will prove (among other things) that this quantity is finite. 

3.2 Bounds of ||p#||^i(z) independent of H 

In this section, we show that n a = supjj e / ,r) k h ^ s finite for a large class of filters, including 
the collection of dilated filters (a m ) m >i of a filter a that will be used in the next section. Recall 
that a m is the filter of length ml + 1 with same order p as a and defined for i = 0, . . . , ml by 

cij/ m if i/m is an integer , 
otherwise. 

As a typical example, if a := a 1 = {1, —2, 1}, then a 2 := {1, 0, —2, 0, 1}. 

Since 7rjy(0) 7^ 0, observe that, for a fixed i 6 Z, the functions iJ 7r#(«) and _ff /0#(«) 
are continuous respectively on [0, 1] and on (0, 1). Moreover, since for any filter a, 

1 lit 

q,r=0,q^r q,r=0 q=0 q=0 

the function H >-» p^j(i) is continuous in 0. In particular, this ensures that for p = 1, ||p a ||£i(z) 
is continuous on [0,1/2). Actually, this may be not continuous in 1/2 but nevertheless n a = 

2 sup ffe [0,1/2] IIpItIU 1 (Z) < +°° for instance /t^ 1,1 ' = 4 and rc't -1 ' 1 ' = 8. We refer to Appendix 
[A] for the computation of the exact values and to Table [3] for the estimation of some other similar 
constants. 

For any filter of order p > 2, observe that 7r"(i) = for all i. Let us consider the following 
assumption on the filter a, denoted H a : 

i 

r a := £ a q a r (q - rf log(|g - r\) ^ 0, (16) 

q,r— 

with the convention 01og(0) = 0. Tab. [2] below shows that Assumption H a is satisfied for a large 
class of filters. Then, from the rule of l'Hospital, 

a ,., E!r=O a <?°'-0?- r + *) 2lo g(k- r ' + *l) 

hm_ p H (i) — j < +°°- 

E g ,r=o a ? a '-(9- r ) 21o g(l9- r l) 

Therefore, under H a , p%(i) is a continuous function of H G [0, 1]. Actually, the same is true for 
the ^-norm of a filter of order p > 2 as stated in Proposition [5] below. 
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a 


1 


2 


3 


4 


5 


P = 


2 


Increments 2 


5.55 


22.18 


49.91 


88.72 


138.63 






Daublets 4 


0.62 


2.47 


5.56 


9.89 


15.45 






Coiflets 6 


0.61 


2.42 


5.45 


9.69 


15.15 


P = 


3 


Increments 3 


13.50 


53.98 


121.46 


215.94 


337.40 






Daublets 6 


0.49 


1.98 


4.45 


7.90 


12.35 


P = 


4 


Increments 4 


41.43 


165.70 


372.84 


662.82 


1035.66 






Daublets 8 


0.45 


1.81 


4.08 


7.25 


11.32 






Symmlets 8 


0.45 


1.81 


4.08 


7.25 


11.32 






Coiflets 12 


0.45 


1.79 


4.03 


7.16 


11.19 



Table 2: Computations of T a ™ for different filters a and its dilatation a' m for to = 1, ... , 5. 

Proposition 5 Let a be a filter of order p > 2 satisfying H a in (|16[) . Then |||0fl-||^i(z) * s o, 
continuous function of H 6 [0, 1] . 

Proof: From ([TT]) . we have 

\J\ 



2H I 2H 



= —i — „ n ^ E a ^ 

q,r— 



For |j| > ^ + 1, we have q — r + j > for < g, r < so that: 

2H 



= ?3 „ „ I |2« 

(?,r= 



b1 2ff 


E £ 


=0 a,ja r |9-r| 2H 








=0 a 9 a r |9-r| 2ff 







1 + J — 

] 

^ (2H)(2H - 1) . . . (2if - fc + 1) fq-r xk 



fc! 

g,r=0 fc=0 



* E — " E <w(<?-r) fe .(r7) 



E„,r=0 a 9 a ''l9- 7 '| 2ff fe=2 P fc!j 



fc=2p J q,r=Q 



Observe that in (IT71) . the outer sum starts at k = 2p. This is due to the property (ITU1) of the 
filter a of order p which implies the following remark: 



k—i 



1 Ik 

^a q a r (q-r) k = ^ a q a r ^ f .j 

g,r— q,r— z— 

- E(-D fc - 4 C)fE^E^H 

= 0iffc<2p-l. 

As a consequence, for p > 2, each summand in the outer sum ([T7|) contains the factor 2if — 2 
in the product (2H)(2H - 1) . . . (2iJ - k + 1). Observe that under H a in (JTBJl, the rule of 
l'Hospital ensures that the function a (H) = (2 — 2H)/(J2 q ^ r Q> q a r \q — r| 2H ) is bounded at 
H = 1~ . Since moreover this function is continuous in we derive, under H a , that ||# a ||oc ■= 
su Pffe[o,i) K(H)\ < 
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Now, from (fTT)) . we have 



< \UH)\\jr-^ %\ k ;fJ f K|K||?-rr* 

fc=0 V y hJ q,r=0 

+00 / I 1 \ k 

\q-r\ 



< \\0 a \Uj\ 2H - 2p £ l« 9 IKH9-H 2p E7fcTT 



(fc + i) v e + i 

q,r=0 k=0 y ' V 

< C(a)\j\ 2H - 2p (18) 
where 



C(a) = ||g.l|ooX:i a »lKllg-r|»f (/ + 1) y / + 1) )< 



Q,r— 



+00. 



When p > 2, the bound (jTSJ) ensures that the convergence of the series X^iez |Pzr(*)l * s un if° rnl 
in H <E [0, 1] and thus H h-> ||^||fi(z) is continuous on [0, 1]. □ 

Proposition [5] proves the following bound is finite for a filter a of order p > 2 satisfying H a : 

k° =2 sup n a H = 2 sup UphIUmz) = 2 SU P IIPhII^MZ) < +°°- ( 19 ) 
He(o,i) He(o,i) He[o,i] 

As a consequence of this result, this means that the constant K a can be obtained by optimizing 
the function H h->- ||fi(z) on the interval [0, 1]. See Tab. [3] below for the computation of such 
constants for different typical filters. 

For dilated increment-type filters, we manage to compute the exact value of ||p#||^i(z) (see 
Appendix [A] for more details) 

fe=l Z^j=l "jJ 2-»j=l 

where ctj = ^ £ ? . r= o a g a r , e(2H — 1) := sign(2H— 1) and where = X)j=o ■ For ^ ne dilated 

g—r=j 

double increments filter a = {1,— 2,l} m for example, this leads to k^ 1 ' -2,1 ^ = 2 x 8/3 = 16/3 



and K^" 2 -!} 2 = 2 x (2 + 25l °i% 2 l los(3) 



7.813554. 
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a 


I 


2 


Q 

O 


4 


u 


P = 


1 


Increments 1 


2 


4 


6 


8 


10 


P = 


2 


Increments 2 


2.667 


3.907 


5.745 


7.565 


9.376 






Daublets 4 


2.250 


4.356 


6.641 


8.906 


11.162 






Coiflets 6 


2.259 


4.327 


6.582 


8.816 


11.042 


P = 


3 


Increments 3 


3.200 


3.783 


5.396 


7.406 


9.200 






Daublets 6 


2.429 


4.516 


6.688 


8.833 


10.966 


P = 


4 


Increments 4 


3.657 


4.304 


6.364 


8.514 


10.350 






Daublets 8 


2.648 


5.026 


7.349 


9.648 


12.044 






Coiflets 12 


2.701 


5.112 


7.459 


9.775 


12.229 



Table 3: Computation of sup ffgJ Wp^P Wt 1 f° r different filters a and for m = 1, . . . , 5. Note that 
I = [0, 0.5] for p = 1 and I = [0, 1] for p > 1. 



4 Confidence intervals of the Hurst parameter 



For any a E (0, 1), denote by q% n (a) :— ((p,, n ) 1 (a;K a ) for • = l,r. In order to make easier the 
presentation, define also 

V/W" 1 / z and Vn a ) :=1+ / 

Note that Remark [2] above ensures that for any a E (0, 1) and for all n > £, x° n _ £ (a) > 0. For 
further reference, observe that for • = I, r and n — > +oo: 



Q?,n-M) - 9 Q («) := V2« a log(l/a). (20) 

In the sequel, we restrict ourselves, to filters of order p > 2 which allows us to make no assumption 
on H . Taking a filter of order p = 1 would have constrained us to assume that if < 1/2. 

4.1 Scaling parameter C known 

In this section, we assume, without loss of generality, that C = 1. Our confidence interval in 
Proposition [6] below is expressed in terms of the reciprocal function of g n (x) :— 2xlog(n) — 
log (7r"(0)), x E (0, 1). In order to ensure that g n is indeed invertible, we assume that 



£q,r=0 a 1 a r lo g(lg ~ r l) |g ~ r ? X ' 



n > exp sup q ' r e ■ . (21) 



In this case, the function g n is a strictly increasing bijection from (0, 1) to ( — log(7TQ (0)), +oo). 
Moreover recall that a filter of length 1 4- 1 requires a sample size n > £ + 1. Obviously, condition 
([2~Tj) only makes sense if the filter a satisfies: 

Eir=o a i a r lo s(k - H)k - A 2x 

sup 1 < +oo. 

xe(o,l) Eor=o a ? a ''l'5 , - r l 2:E 
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a 


1 


2 


3 


4 


5 


P 


= 2 


Increments 2 


3 


4 


6 


9 


11 






Daublets 4 


4 


6 


10 


13 


15 






Coinets 6 


6 


11 


15 


21 


2G 


P 


= 3 


Increments 3 


4 


6 


10 


13 


15 






Daublets 6 


6 


11 


15 


21 


2G 


P 


= 4 


Increments 4 


4 


9 


13 


17 


21 






Daublets 8 


7 


15 


22 


29 


3G 






Symmlets 8 


7 


15 


22 


29 


3G 






Coinets 12 


12 


23 


34 


44 


56 



Table 4: Minimal sample size n required to satisfy (|2ip for different dilated niters a m of different 
orders p. 



Since \iuy x ^i- r=0 a q a r \q — r\ = 0~ (we stress that this function vanishes with non-positive 
values of because it is continuous, negative in x = 0, see (|15[) . and does not vanish), the previous 
condition is equivalent to the more explicit following one 



a q a r \og(\q - r\)(q - rf > 0. 

q,r—0 



(22) 



Table 0] exhibits the minimal sample size n required to satisfy ([2~Tj) for different filters a rn (for 
m = 1, . . . , 5) with different order p = 2,3, 4. Obviously, condition (|2"2"j) is in force for all these 
filters. 

We state now our main result when the scaling parameter is known: 
Proposition 6 Let a € (0, 1) be fixed and a be a filter satisfying H a in (|16l) 

1. For n > £ + 1, we have: 

P( log (xt >n _t{a/2)) - log < g n (H) < log « n _,(a/2)) - log ) > 1 - a. (23) 

2. Moreover if the filter a satisfies (|22[) and n > £ + 1 satisfies (|21|) , we have: 



whe 



He 



max 



> 1 - a, 



(24) 



(o^-^log (xl n _ e (a/2))~log(S a n )) 



H^(a) := min {r,g~ l (log «„_,(a/2)) - log (££)) 

5. is n-> -f-oOj i/ie proposed confidence interval in (|24p satisfies almost surely 

'Hf(a),H s n ^(a)]^{H} 
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and the length fi n of the confidence interval satisfies 

2q a {a/2) 1 q a (a/2) 



fin 



where q a is defined above in (|2"CT 



Vn 9n{ H ) \/nlog(n)' 



Remark 3 PronositionXfh generalizes Provosition\l\ derived from Breton et al. 1 200& ). The scal- 
ing parameter is still assumed to be known. However, we do not need to know an upper-bound 
of H and our condition on n is much sharper than the one required in Proposition [71 As an 
example, for a — (1, —2, 1), condition (|21[) is satisfied for all n > 3, whereas the minimal sample 
size allowing to derive a confidence interval from Proposition[l\is 1108 for a = 5% and H* = 0.8. 

Proof: Consider the set 

A := {-ql n _ e (a/2) < yfc=£V° < <£ n _<(a/2)} . 

The bound (IT3|) entails P(A) > 1 — ^ — ^ = 1 — a - It is now sufficient to notice that 



A = K n _,(a/2)<l + ^<< n _,(a/2)} 



n 2H 



l0 § On ^ 9n(H) < log 



Ca I ° I Co 

/ \ u n 

which proves ((23)) . Next, since under (|2"Tj) and (|2"2")l . g n is an increasing bijection, (|2~31) comes 
immediately from (|2"3")l . Finally, from (I2U1) . we have 

log(xl n _ e (a/2))~-^B and log « n _,(«/2)) ~ Q -^B 
v n V 71 

asti^ +oo. Moreover, since 1 + V£ = ^m S n = S^e 9n(H \ using the LLN in (12]), we have 
almost surely 

- log (S£) = - log (1 + V£) + g„(H) = g n (H) - V£(l + o(l)) ~ 

It is proved in ICoeuriollvl (|200lh (Proposition 1) that V£ converges almost surely towards for 
any filter and for all H £ (0, 1) which implies the almost sure convergence of the confidence 
interval and the asymptotic behavior of the length fi n of the confidence interval. □ 



4.2 Scaling parameter C unknown 

The idea to construct confidence intervals when the scaling coefficient C is unknown consists in 
using the collection of the dilated filters a m defined in (|T3|) . 

Let us first introduce some specific notation: let M > 2 and consider a vector d = (di, . . . , dM) T 
with non zero real components such that Yli=i di = and such that cPLm > 0, where 
Lm = (log(m)) =1 M . Denote by I~ and I + the subsets of {1, . . . , M} defined by 

I~ = {i £ {1,...,M} : di < 0} and 1+ = {i £ {1, . . . ,M} : di > 0} . 
The following confidence interval is expressed in terms of Ls„ := (log (S 1 " )) _i m- 
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Proposition 7 Let a G (0, 1) be fixed and denote by L X inf and L x »up the two following vectors 
with components 



log («/2M)) ifmer 



(Lvu) = 

log [x a r 2_ ml (a/2M)J i/m e J+ 
i. Let n > M£ + 1. T/ien we fccrae 



> ( L x= up ) m - 



log fc-^ («/2M)) i/mer 
log K™ « (o/2M) ) i/me/+. 



> 1 - a 



(25) 



where 



H™ f (a) = max |^0, 
^ up (a) = minfl, 



1 



2cFLm 
1 



(d T U 



d T L 



;d T L s „ - d T L x - P ) 



2cFLm 

2. As n — > +oo ; i/ie proposed confidence interval in (|25[) satisfies almost surely 



and its length /i„ satisfies 



(f (L X i°f -L x .up) 1 d T q M (a/2M) 



where qjvi(a/2A/) is t/ie vector of length M with components defined by 



(q M (a/2M)) m := 



u>i£/i g a defined in (|2CT 



-g am (a/2Af) if me I~ 
q am (a/2M) if m € 7+ 



Remark 4 Proposition generalizes Proposition since t/iis new confidence interval does not 
assume that the scaling parameter, C is known. More specifically, note that the definition of the 
interval does not depend on C. Note also, that if Bh were not observed on [0, 1) but with a 
dilatation factor, then the confidence interval would remain unchanged. 

Proof: For m = 1, . . . , M, we consider the following event 

A m := {aC-mi («/ 2M ) < 1 + K m < KZ-mi («/2M)} . 
The bounds (O entails that P(A m ) >l-jt~jt=l-jj- First, recall that 



ya 



C%^(0) 



S a n -l= 7 x 



1 



,5° - 1 with 7 := 7c,ff,n = 



C^(0)' 



The crucial point in the definition of the confidence interval relies on the fact that 7 is independent 
of m. Second, note that for m = 1, . . . , M: 
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< log 1 



Kf) <log(<:_ m ,(a/2A/))} 



= {^g (x?Z-mt («/2M) 
= {log ^;_ mJ ( a /2M)) - log( 7 ) < log (sf) - 2H log(m) 

< log («:_ m ,(a/2Af)) -log( 7 )} 
= {log (s ? <f ) - log (<:_ mf ( a /2M)) + log( 7 ) < 2# log(m) 

< log (Sf) - log ((xf n _ mt {a/2M)) + log( 7 )} 

= {d m ((L s J m - (L x ,..) m + log( 7 )) < 2d m H(L M )m < d m ((L s J m - (L x -p) m + log( 7 )) } . 

Next, we consider the following event 

B := {d T L Sn - d T L x »f + d T l log( 7 ) < 2Hd T L M < d T L s „ - d T L x - P + d T l log( 7 )} 
= {d T L s „ - d T L X mf < 2ffd T L M < d T L s „ - d T L x - P } 

where 1 = (1, . . . , 1) T . Since A x n A 2 n . . . n A M C 5, setting A c = \ A, we have 

P(B) > p(Ai n . . . n A M ) = i - p ((Ai n . . . n A M ) C ) = l - v(A$ u . . . u A C M ) 

M M 

> l-£p(^) = £P(4»)-(M-1) (26) 



-oo, the following estimate 



> M(l-— J-(M-l) = l-a, 



which ends the proof of 
holds almost surely 



Next with the LLN in fl2"]h as 



log(sf) = 2 J fflog(m)-log( 7 )+log(l + V; Q 

= 2 J fflog(m)-log( 7 ) + K m (l + (l)) 
and implies that almost surely, when n — > +oo, 



d T L Sn 



2#d T L M - d T l log( 7 ) + d T (V£ m ) (1 + o(l)) 

\ / m=l,...,M 



2H d T L M + d T 



ya 



m=l,...,M 



(l + o(l)) 



2Hd T L M - 

From (|2"U)) . one has also the following estimates as n — > +oo: 



~q am {a/2M) if m G J~ 
q a (a/2M) if m € /+ 



(L x s«p) ~-=: 



q am (a/2M) if m G J" 



\ -<f (a/2M) if me/ 4 



These different results imply the almost sure convergence of the confidence interval towards 
{iJ}. For the asymptotic of the length \x n of the confidence interval, it is sufficient to note that 
(L xi „ f - L XB u P ) ~ -^q M (a/2Af). □ 
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5 Simulations and discussion 



5.1 Confidence intervals based on the central limit theorem 
5.1.1 Methodology 



There exists a very wide litterature on the estimation of the Hurst parameter, see e.g. ICoeuriollv 



(2000) and references therein. For all of the available procedures, the confidence interval comes 



from a limit theorem so that it is of asymptotic very nature. In contrast, our confidence intervals 
in and (1231) are non-asymptotic since they are based on concentration inequalities. In order 
to compare our procedures, we choose to focus only on one of these procedures which has several 
similarities with this paper. These procedures are based on discrete filtering and are presented 



in detail in ICoeuriollvl (|2001l) . For the sake of sclf-containess, we first summarize them: 

• Scaling parameter C known. The procedure is based on the fact that almost surely 
tt" (o) ffit ~~ * 1' n +°°- With the same function g n {x) — 2x log(n) — log(7r£(0)) as the one 
used to derive the confidence interval in Proposition [6j this yields the estimator: 

K td (a) -^(-logOS*)). 



Note that the confidence interval (|24| is very close to this estimator. In particular, the 
middle of the interval ([24]) behaves asymptotically as H^ td (a). 

Scaling parameter C unknown. The idea of ICoeuriollvl (<2000h in this context is to 



use the following property of quadratic variations of dilated filters E[S^ ] = m 2H ^f with 

7 : = — n 2H and the almost sure convergence of S% /E[S^ ] towards 1 for all m. The 
idea is then to estimate H via a simple linear regression of Ls n on 2Lm for M dilated 
filters. Here, the notation Ls„ and Lm are the same as the ones in Proposition [7] This 
leads to the estimator 

Hr(a,M) :=|pqp 
where A = (log(m) — J2m=i l°g( m ) ) ■ There is again an analogy between this 

V / m—l....,M 

estimator and our confidence interval in Proposition [7J Indeed, with d = A, the interval 
in (|2"5j) rewrites 

' A T (L Sn - Lxjjf) \ . / A T (L S „ - L X . UP 
0, T^r-r^ — - — ,mm 1, 



2||A|| 2 i' I ' 2||A|| 

since d T Ljvi = A T A = ||Aj| 2 . Again, the middle of this interval behaves asymptotically 
as H^ en (a, M). In the particular case M = 2 the estimator H^ en (a : 2) takes the simple 
following form 

and the bounds of the interval in (1251) rewrite as 





1 'l 0g f ^ i- ( <n-2i(a/±) 



2 log 2 \ ° °\xf n _ e (a/4) 




1 l log l^ \ ,_J*?, n -it(<*/4) 



2 log 2 \ I 5« M °l (a/4) 




1G 



5.1.2 Asymptotic confidence intervals 

We refer the reader to Coeuriollv (2001) where the following central limit theorems (CLT) are 
proved for H^ d (a) and H% en (a,M) 



/ nlog(n) 



(27) 



where — ► stands for the convergence in distribution, Af(0, 1) is the normal standard distribution 
and ° 2 s td( H ) ■= 3 HPhII^(z)> and 



hi— — ^ ► JVYO.l), n^+oo 

(T s td{Hii ,M) 



(28) 



where <Tg en (H, M) := -jpqjT- where G is the (MxM)-matrix defined by G mii 
for mi, ?72 2 = 1, . . . , M, and for all i €E Z 



2(Z) 



a 1 ,a 2 / -\ 



VfllT + i 



2H 



Note that in the special case where M — 2, the constant al en (H, 2) takes the simple form 



2(log2)2 



Ph 



£ 2 (Z) 



Ph 



£ 2 (Z) 



Thanks to the CLTs, (|2"T|) and (j2"5)l an asymptotic confidence interval to the level 1 — a, a € (0, 1), 
can be easily constructed 



ICf(a) 



a' 



0,H--$-\l-a/2) x — ,min [l,H' + <p-\l- a/2) x — 



(29) 



where • = std,gen, v^ d = y/n\og(n), vff n = y/n and <E> is the cumulative distribution function 
of a standard Gaussian random variable. 



5.2 Comparisons of approaches 

In the following tables, we compare, via Monte-Carlo experiments, the confidence intervals based 
on concentration inequalities (|24l) , (|2"5|) and on central limit theorems (f2U|) . The fractional Brow- 
nian motions have been generated by using the circulant matrix method (e.g. Kent and Wood! 



(|l997h . ICoeuriollvl tQ0(t) ). We have realized a very large simulation study. The "best" results 



(in terms of choices of the filters a, of the maximum dilatation factor M) are summarized in 
Table [S] for the standard fractional Brownian motion (i.e. C = 1) and in Table [5] for the general 
one (i.e. C unknown). 

In Figure[Tl we also compare, in terms of H, the asymptotic lengths of the confidence intervals 
obtained by each approach. 
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H = 0.2 






H = 0.5 






H = 0.8 








Cover. 


Length 


H 


Cover. 


Length 


H 


Cover. 


Length 


H 


n = 50 


CI[i2] 


100.0 


0.2191 


0.1875 


100.0 


2029 


0.4832 


100.0 


0.1553 


0.7824 




CLT[i2] 


95.2 


0.1330 


2058 


97.0 


0.1227 


0.5013 


99.6 


0.1125 


8003 




CI[d4] 


100.0 


2086 


0.1886 


100.0 


0.1941 


0.4841 


100.0 


0.1482 


0.7834 




CLT[d4] 


94.6 


0.1217 


0.2050 


97.2 


0.1133 


0.5004 


99.2 


0.1076 


0.7999 


n = 100 


CI[i2] 


100.0 


0.1298 


0.1936 


100.0 


0.1212 


0.4946 


100.0 


0952 


0.7931 




CLT i2 


95.0 


0800 


2009 


97.6 


0.0737 


0.5017 


99.8 


0.0676 


8003 




CI[d4 


100.0 


0.1224 


0.1941 


100.0 


0.1149 


0.4949 


100.0 


0902 


0.7933 




CLT[d4] 


95.6 


0.0732 


0.2005 


96.4 


0.0680 


0.5012 


99.6 


0.0646 


0.7997 


n = 500 


CI[i2] 


99.6 


0.0430 


0.1994 


100.0 


0.0408 


0.4988 


99.8 


0336 


0.7988 




CLT i2 


94.4 


0265 


0.2004 


96.4 


0.0244 


0.4998 


98.8 


0.0224 


0.7998 




CI[d4 


99.6 


0.0402 


0.1995 


100.0 


0383 


0.4990 


99.8 


0.0316 


0.7989 




CLT[d4l 

L J 


95.4 


0.0243 


0.2003 


96.0 


0.0225 


0.4999 


98.4 


0.0214 


0.7998 


n = 1000 


CI[i2] 


100.0 


0.0274 


0.1998 


100.0 


0.0262 


0.4996 


100.0 


0.0219 


0.7997 




CLT[i2] 


96.6 


0.0169 


0.2003 


97.6 


0.0155 


0.5000 


99.2 


0.0142 


0.8001 




CI[d4] 


100.0 


0.0256 


0.1998 


100.0 


0.0245 


0.4996 


100.0 


0.0205 


0.7998 




CLT[d4] 


96.4 


0.0154 


0.2002 


97.2 


0.0143 


0.5000 


98.8 


0.0136 


0.8001 


n = 10000 


CI[i2] 


99.8 


0.0066 


0.2000 


100.0 


0.0063 


0.4999 


100.0 


0.0055 


0.8000 




CLT[i2] 


94.2 


0.0040 


0.2000 


96.2 


0.0037 


0.5000 


98.4 


0.0034 


0.8000 




CI[d4] 


99.8 


0.0061 


0.2000 


99.8 


0.0059 


0.5000 


100.0 


0.0051 


0.8000 




CLT[d4] 


94.4 


0.0037 


0.2000 


95.0 


0.0034 


0.5000 


98.2 


0.0032 


0.8000 



Table 5: Monte-carlo experiments based on 500 replications of a fractional Brownian motion 
with Hurst parameter H = 0.2, 0.5, 0.8 and scaling coefficient C = 1 (assumed to be known) and 
for different values of the sample size n. The filters i2 and d4 denote respectively the filter of 
Increments of order 2 and the Daublets 4. 
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Table 6: Monte-carlo experiments based on 500 replications of a fractional Brownian motion 
with Hurst parameter H = 0.2, 0.5, 0.8 and scaling coefficient C = 1 (assumed to be unknown), 
for M = 2, 5 and for different values of the sample size. The filters i2 and d4 denote respectively 
the filter of Increments of order 2 and the Daublets 4. For these simulations the vector d has 
been fixed to the vector A. 
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5.3 Discussion 



We propose non-asymptotic confidence intervals for the Hurst parameter of a standard or non- 
standard fBm based on concentration inequalities. They are computable in particular for small 
sample size and several theoretical improvements are obtained: 

• W hen the scaling para meter C is known, we have refined the confidence interval proposed 



m 



Breton et al. (2009): the upper bound H < H* < 1 is relaxed, the condition on the 



sample size n is sharper and our new confidence intervals are valid for a large class of 
filter a. 

As a by-product in our way to optimize t he nu meric bounds, we have slightly improved the 
bounds obtained bv lNourdin and Viens (|2009l) in the general concentration inequality (see 



Proposition [2]). 

The case where C is unknown has never been considered with concentration inequalities 
before Proposition [7] 

• The asymptotic properties are similar to that of confidence intervals based on central limit 
theorems. More specifically, the length of the confidence intervals derived by concentration 
inequalities behaves asymptotically as the ones of confidence intervals based on central 
limit theorems, that is l/(y / nlog(n)) when C is known and ^/\/n when C is unknown. 

The comparison with confidence interval based on CLT is contrasted: while the Monte- 
Carlo experiments are correct when C is known (in terms of coverage rate and of lengths of the 
confidence intervals), they are not good when C is unknown: the lengths equal often 1, i.e. the 
intervals correspond to (0, 1), when the sample size is small and are about five times larger when 
n is large. In fact, the confidence intervals derived from concentration inequalities are too much 
"sympathetic": the coverage rate is rather far from 1 — a (based on 500 replications, it is even 
often equal to 100%). From a statistical point of view, this is the main reason why the length of 
the confidence interval is sometimes much larger than the ones based on central limit theorems. 
From a mathematical point of view, this is due to the fact that, in Proposition!?! the dilatations of 
a filter are actually handled separately. As a consequence, the errors induced by each dilatation, 
and controled by the concentration inequalities ©-Q, add up, see (|2l))) . This explains that the 
proposed confidence interval based on concentration inequalities are less performing in this case 
while, in comparison, multivariate CLT are used for standard confidence intervals. Improvements 
would require to use multivariate concentration inequalities, generalizing Proposition [2J which, 
at the moment, are not available. This is the aim of future research to obtain such improvements. 

As a conclusion, this work is the first attempt to define computable confidence intervals for 
the Hurst parameter H of a standard and a non-standard fractional Brownian motion with an- 
ot her approach than the classical one based on central limit theorems (at the very exception 
of Breton et al. where the first non- asymptotic confidence intervals were derived for the 



standard fBM with a more theoretical motivation). We did not get around the question of the 
numerical performances via Monte-Carlo experiments. The conclusion is that, based on concen- 
tration inequalities, confidence intervals can be proposed for a large class of filters and without 
assumption on the Hurst parameter. The performances are comparable to the stantard confidence 
interval based on CLT when the scale parameter C is known, while the procedure is underper- 
forming when C is unknown. This later case requires preliminary theoretical improvements for 
multivariate Gaussian quadratic forms that motivate our future studies. 
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Ratio of asymp. lengths for Increments 2 



Ratio of asymp. lengths for Daublets 4 




Ratio of asymp. lengths for Increments 2 Ratio of asymp. lengths for Daublets 4 




0.0 0.2 0.4 0.6 0,8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 1: Ratio of asymptotic lengths of confidence intervals of procedures derived by concen- 
tration inequalities and central limit theorem when the scaling parameter C is known (top) and 
unknown (bottom). The confidence level equals 1 — a = 95%. For the general procedure, the 
vector d has been fixed to d := Lm — Lm 
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A Exact computations of 4-norm for filtered fBm 

In this section, we describe how explicit exact bound can be obtained for the correlation of a 
filtered fBm. Let a be a filter of order p and length I. Its covariance function is given by 

t t 
= ~\ E a q a r \q-r + k\ 2H = £ a,\j + k\ 2H 

q,r=0 j=—i 

where otj = Y^q,r=o a q a r . Note that 

q—r=j 

• aj — at-j, in particular 7r|j-(0) — — X^=i a jj 2H \ 

• T,j=-i a 3 = Yfq,r a q a r = 0, 



• for all h < 2p — 1 , we have 



E i h(x i = E i h E a i ar = E E {z-r) h a q a r = ^ (q - r) h a q a r 
j=-e j=-e q-r=j j=-£q-r=j q,r=o 



k 

q,r=0k=0 v 



0. (30) 



X^o a i = ~ a o = - J2 q =o a l < 0, on = a a e . 
A crucial observation is that, at least for \k\ large enough, all the irfj(k), and thus all the p a (k), 
have the same sign. Indeed, using (|3U)) . we have for |fe| > I: 



n%(k) = -i^a,(|fc +i r + |fc-ir-2|fcn 



3 = 1 

j=i 

:,, 2H ^ I (2ff)(2g-l)...(2g-2i + l; , v 

= y pp* 12." 

(2ff)(2ff-l)...(2ff-2p+l) | ^ 2p 

This observation allows to reduce the computation of the ^-norm |||0|f||^i(z), which is an infinite 
sum with modulus, to an infinite sum of correlations but without modulus plus some finite sum 
(with modulus remaining). Essentially, it remains to compute the sum of correlation without 
modulus. This is done below. But observe first that if there exists some k(H, a) £ N so that the 
correlations p%(k) have all the same sign for \k\ > k(H,a) large enough. The value k(H,a) is 
not known in general. However for some family of filters (including increment- type filters in and 
their dilatations (in) m , n, m > 1), k(H,a) is known and explicit computations are tractable: 
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Proposition 8 For a dilated increment-type filter a £ {(in) m : n,m> 1}, we have k(H,a) = £, 
i.e. the following property holds true: 

for all \j\ > £, -K a H {j) is of the same sign as (-1) P+1 (2H - 1). (31) 

Proof: Let 9 m (f)(x) = f(x + m) — 2f(x) + f(x — rn). Observe that if / is a convex 
(resp. concave) function, then 9 m (f)(x) > (resp. O m (f)(x) < 0). For the il hlter, we have 
n l ^(x) — ^9i(\x\ 2H ), for the i2 filter, we have tt 1 h{x) — — ^9° 2 {\x\ 2H ) and more generally for the 

m-dilatation of the in filter, we have tt^^ (x) — — % — 9™(\x\ 2H ). 

Observe also that the function \x\ 2H and all its iterated derivatives (\x\ 2H Y 2p ^ of even order 
are convex if if > 1/2, concave if H < 1/2. By an immediate induction on n, we show that the 
same holds true for all 9™(\x\ 2H ). In particular for \j\ > £m, we obtain that t^™' 1 (j) is of the 
same sign as (-l) n+1 (2H - 1). □ 

Obviously, the property ([3~Tj) does not hold true for any filter (consider for instance {1, —4, 5, —2}). 
In order to make easier our following explicit computation to derive exact value for ||p a ||^i(z)j we 
consider a filter a satisfying (|3"Tj) but we stress that for each particular filter the same strategy 
applies with some specific k(H,a). First, for all N > I, we have: 

N N £ 

j=e j=i k=—i 

£ N £ N+k 

k=-l j=l k=- 
£ 

= a~iS§_ e + E a k (S§ +k - Sf +h _x t 

k=-l+l 

£ I N+k 

= a. t S§_ t + E <*k[S§-i+ E \i\ 2H -Sf +k -x 
k =-e+i \ j=N—t+i 

C£ \ ( 1 N+k \ / e x 

E«* ) s %-i+[ E «* E i;r - E «^Vi 
k=-e / \k=-e+i j=N-e+i J Vfc=-^+i / 

£ 

= X N - E a kS" +k _i 
k=—i+l 

where = Ylj=a 3 2H an ^ 

N+£ 

XN = 




«\?h I ( 2Hi 2H(2H-l)i 2 / 1 \\ A \ , s 
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But 



and 



E E afc = E {t + k)a k = (£ + k)a k = 



i=0 k=l-i 



21-1 



E *E 

j=0 \ fe=i— . 



E 



/ 2*-l \ t I 21-1 

a >* E = E afe E 



=i-As 



k=-e 



e-k 



E a * ( 



/2£(2^-l) fc- 1) 



k=-e 



because of l[30]). We obtain zjv = 0((iV + f) 2ff ~ 2 ) -> 0, iV -> +oo. Actually, expanding 
(1 — i/(iV + ^)) 2ff to the (2p — l)-th order in (j3"2"j) . and since Ei=i i s a polynomial in N of 
degree + ([301) shows that x n = 0((N + £) 2H - 2p+1 ) . Finally with the property (JSTJ), we have: 

+oo I 

2Ei^wi = (- 1 ) p+le (( 2i? - 1 ) E «^Vi 

and 



-t-oo 



k=l 



k=t 



/c=i 



2ff 



2 if 



i+E- 



2 if 



fc=l 



2/f 



(-lf +1 e(2ff- 1) 



(-lf +1 e(2ff- 1) 



2ff I 



El qH 
k=-i+l akD t+k-l 



£,=1^.? 



2ff 



(33) 



where we recall that e(2H — 1) = sign(2H — 1). First, note that the modulus has been removed 
in the denominator of ((33)) according to the following observation: 



e i 

e^ 2H 4e 



j=l 



o=-t 



1 

2 E a J : = E a J ~ a ° 

j=-< 



-a < 0. 



•2_ff 



Since we assume moreover 7r^(0) 7^ 0, this means that 7T^-(0) > and that J2j=i a j3 

Next, note that (|33p is an explicit expression involving only finite sums and can be easily 
explicitely optimized for H 6 (0, 1) for every given a satisfying H a . Note that, for p > 2, when 
H — > 1, right-hand side of ((33]) remains well defined. Observe first that since for any fixed k, 



lim ff _,i Sf = Sl = fc(fc+1) 6 (2fc " 1) , we have using Qty 



£ 1 * 

hm ^ QfcSj* fc _ x = - E afcO? + fc-l)(^ + /c)(2£ + 2£:-l)=0. 



fe=-£+l 
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The same holds true for Y^,j=-i a j\j + k\ 2H an d 2i=i a jJ 2H > but under H° in (H"6"l) . the rule of 
l'Hospital entails \im H ^i- fl/off ||.£i(z) exists and is finite. Since obviously, is a continu- 

ous function of H £ [0, 1), this ensures the continuity of ||pf;-||,«i(z) on [0, 1] and the constant K a 
in our confidence interval is obtained by maximazing the explicit function in (|33[) . 

Dilated simple increments (il) m = {—1,1}™. In this case, t = m, p = 1, atj = for 

1 < j < m and olq — 2, a± m — — 1 so that (|33[) rewrites: 



{-MF 
Ph 



m—1 



|j + m| 



2/i 



2|j r W + |j-m| 2ff | 



oH 

°2m-l 



ZJ m-l 



For instance for m = 1, 



3=1 

I I- 1 - 1 }! 



«i(Z) 



-,2/f 



2 and for m 



,2ff 



(34) 



Ph 



2 4"+9 f ff -l i SQ that 



K 41 = 4 and = 8 (recall that in this case, we optimize for H 6 (0, 1/2]). 

In general, since the right-hand side of (f3~4"|) is a continuous function of H, and since for all 

= 2m while 



k > 1, Si 



1/2 _ fc(fc+l) 



we have lim 



ff-Ki/2)- 



{-i,i}" 
Ph 



£!(Z) 



Pl/2 



exhibiting a discontinuity of the ^ 1 -norm for the dilated i 1 filters. 

Dilated double increments (i2) m = {1,— 2, l} m . In this case, £ = 2m, p = 2 and ao = 6, 

a± m = —4, a±2m = 1, «j = 0, j 7^ 0, ±m,±2m, so that ([33)) rewrites: 



{1,-2,1}™ I 
IPjf IUl(Z) 



2m-l 

+ E- 

fe=i 

-e(l-2ifl 



|/s-2m| 2H -4|fc- 



|2H 



6|fc| 



2H 



4I&- 



|2H 



2m 



2ff I 



,2H 



-45 ff 



45. 



rn 
•l 

3m-l 



(4 - 4 H ) 



if 



5 



m 2ff( 4 _ 4 H) 



In order to obtain explicit values, we focus on the cases m = 1 and m — 2. First, for m = 1, (|33[) 
reduces to 



I ^ 2 II 

\Ph\\i i (z) 



10-7x4" +2x9' 
4=4~S 



H < 1/2 
H > 1/2 



and elementary computations entail: 



Next, for m = 2, since 

r (*2) 2 



2^T (i) 

24^(3) 



2tt^ (2) 



expression (1331) reduces to 



2x lim llp'il 

7^0+ 



= 2 1 



5\ _ 16 

3 J ~ y* 



-2 + 3x9 ff - 25 ff 
-7x4 H +4x 16 ff -36 ff 
3 - 6 x 9 H + 4 x 25 H - 49 H 



> ViJ e (0, 1) 
<0 ViJe(0,l) 
<0 V# 6(0,1) 



(i2)' 



I ^ _|_ -B+1UX4" +12XU" -7Xlt> -»X25 +2X3ti +2X49 £ or < J/2 

z) = { 1 + - 4+4x4H \^ljf-^ for if > 1/2. 

An elementary study of this function, together with the rule of l'Hospital, entails that 



2x sup 
-He [o,i] 



{i2)< 

Ph 



^i(Z) 



2 x lim 

fl->4- 



(i2)' 



£!(Z) 



2 1 



251og(5) - 271og(3) 
81og(2) 



7.813554. 
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